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SUMMARY 


In recent years we have developed high accuracy finite 
difference approximations for partial differential equations 
of elliptic type, with particular emphasis on the 
convection- diffusion equation. These approximations are of 
compact type, have a local truncation error of fourth order, 
and allow the use of standard iterative schemes to solve the 
resulting systems of algebraic equations. In this paper , we 
extend these high accuracy approximations to the solution 
of Navier- Stokes equations. Solutions are obtained for the 
model problem of driven cavity and are compared with 
solutions obtained using other approximations and those 
obtained by other authors. It is discovered that the high 
order approximations do indeed produce high accuracy 
solutions and have a potential for use in solving important 
problems of viscous fluid flows. 


*Work funded by Space Act Agreement C-99066-G. 



INTRODUCTION 


The basic model for the fluid dynamics phenomena 
representing two and three dimensional flows of an 
incompressible viscous fluid is given by the Navier- Stokes 
equations that represent the conservation of mass, momentum 
and energy. These equations are highly nonlinear and are 
very difficult to solve, especially when the approximate 
solutions are required to have a high accuracy. A related 
problem is that of obtaining highly accurate solutions of 
the convection- diffusion equation especially when 
convection is the dominating phenomena. 

Prior to the last decade, many researchers examined a 
number of first and second order finite difference and 
finite element methods that were accurate and stable. In the 
area of finite difference methods, it was discovered that 
although central difference approximations were locally 
second order accurate they often suffered from computational 
instability and the resulting solutions exhibited 
nonphysical oscillations. The upwind difference 
approximations were computationally stable though only first 
order accurate and the resulting solutions exhibited the 
effects of artificial viscosity. The second order upwind 
methods were no better and the higher order finite 
difference methods of conventional type were computationally 
inef f ic ient . 
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An exception has been found in the high order finite 
difference schemes of compact type that are computationally 
efficient and stable and yield highly accurate numerical 
solutions at least for the linear and quasilinear partial 
differential equations. Simplest version of such compact 
schemes is given for the Poisson equation 

9 2 u/3x 2 + 9 2 u/9y 2 = f(x,y) 

which can be discretized at a grid point (x,y) by a nine 
point finite difference approximation: 

4[u 2 + u 2 + Uj+ u t ] + u 5 + u 6 + u 7 + u,- 20u o 

= h 2 [ f ! + f 2 + f,+ f,+ 8f o 1/2 . 

(See Figure 1 for the computational stencil). This 
approximation was named Mehrstellenverf ahren by Collatz [3] 
(see also [10]). It has a local truncation error of order 
h 4 and is an approximation of compact type as it involves 
only the eight nearest neighbors of the point (x,y). This 
type of approximations have been obtained for other elliptic 
equations by many researchers: the Hodie schemes of Lynch 

and Rice [14] (see also [2]), the O.C.I. schemes of Berger 
et al.[l], and the SCHOS schemes of Gupta et al. [12,13] all 
reduce to the above difference approximation when applied to 
the Poisson equation. Similar compact schemes of order h‘ 
have also been obtained [15]. 
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The compact schemes of Gupta et al. (called SCHOS in the 
earlier papers) were applied to the convection-diffusion 
equations in particular and were found to yield high 

accuracy when applied to a large number of test problems 
including problems of convection dominated flows [12]. In 
this paper, we extend these finite difference schemes to the 
Navier- Stokes equations. As a test of this method, we solve 
the model problem of lid driven cavity for small to moderate 
values of the Reynolds number and compare our numerical 

solutions with the highly accurate benchmark solutions 

available in the literature. We find that our method yields 
high accuracy even though we use a relatively coarse grid. 

In the next section we describe the fourth order compact 
difference schemes for the convection- diffusion equation 
and for the Navier- Stokes equations. The model problem of 
the lid driven cavity is described in the following section 
that follows with detailed comparisons of our solutions with 
the existing solutions in the literature. 

FINITE DIFFERENCE APPROXIMATIONS 

Consider the following steady convection- diffusion 

equation: 

3 2 u/3x 2 +3 2 u/3y 2 +p(x,y)3u/3x +q(x,y)3u/3y =f(x,y). (1) 
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Finite difference approximations of eq.(l) at the grid point 
(x,y) are obtained in terms of the function values of u, p, 
q, and f at (x,y) and its neighbours. Assuming a uniform 
grid in both x and y directions, we number the grid points 
( x , y ) , ( x+h , y ) , (x,y+h), (x-h,y), (x,y-h), (x+h,y+h), 

(x-h,y+h), (x-h,y-h), (x+h,y-h) as 0 , 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 

respectively (see Figure 1). In writing the finite 

difference approximations a single subscript * j ' denotes the 
corresponding function value at the grid point numbered 'j'. 

The usual central difference approximation of eq.(l) at 
the point (x,y) is given by 

[ u, + Uj + u 3 + u 4 - 4u 0 ] 

+p 0 h( u 3 - u 3 )/2 + q 0 h( u 2 - u 4 /2 = h 2 f 0 . (2) 

This approximation is obtained by replacing all derivatives 
in eq.(l) by central differences that have truncation errors 
of second order. A high accuracy approximation of eq.(l) at 
the grid point (x,y) is given by 

C 1 U 1 + C 2 U 2 +CjU, +c 4 u 4 +c 5 u 5 +c 4 u 4 +c,u, +c.u, -c 0 u 0 
=h 2 [f 2 +f 2 +f 3 +f 4 +8f o ]/2 + h 3 [p 0 (f 1 -f 3 )+q 0 (f 2 -f 4 ) ]/4, (3) 

where , 

c 2 = 4 + h/4(4p 0 + 3p 1 +p 2 -p 3 +p«) 

+h 2 /8 [ 4p 0 2 + p 0 (pi-p,) + q.(Pa~*P«)]' 

c 2 = 4 + h/4(4q 0 +q 1 +3q 2 +q 3 -q 4 ) 

+h 2 /8 [ 4q 0 2 + p 0 (qi-q 3 ) + q.fq^qJ], 
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c, = 4 - h/4(4p 0 -p 1 +p 2 +3p 3 +p 4 ) 

+ h 2 /8 [ 4p 0 2 - p o ( p i ~p j ) - q„(p 2 -p,)], 

c 4 = 4 - h/4(4q 0 +q 1 -q 2 +q 3 +3q 4 ) 

+ h 2 /8 [ 4q 0 2 - p 0 (qrqj - q,(q 2 -q,)], (4) 

c 5 = 1 + h/2(p 0 +q 0 ) + h/8 (q 1 +p 2 -q 3 -p 4 ) + h 2 /4 p 0 q 0 , 

c 4 = 1 - h/2 (p 0 -q 0 ) - h/8 (q ! +p 2 -q 3 -p 4 ) - h 2 /4 p 0 q 0 , 

c 7 = 1 - h/2 (p 0 +q 0 ) + h/8 (q 1 +p 2 -q 3 -p 4 ) + h 2 /4 p„q 0 , 

c 8 = 1 + h/2 (p 0 -q 0 ) - h/8 (q 3 +p 2 -q 3 -p 4 ) - h 2 /4 p 0 q 0 , 

c 0 = 20 + h 2 (p 0 2 + q 0 2 ) + h (p 3 -p 3 ) + h(q 2 -q 4 ). 

This approximation is of compact type as it involves only 
the nine grid values of u at the point (x,y) and its eight 
nearest neighbors (see Figure 1). This approximation has a 
local truncation error of order h 4 . Detailed derivation of 

this approximation is given by Gupta et al. [12], (This 

approximation has recently been rediscovered by Dennis and 
Hudson [21].) Similar high order approximation for general 
second order elliptic equations are given in [13]. Other 

compact approximations of this type have been obtained for 

the Poisson equation [10], the Helmholtz equation [15], and 
the biharmonic equation [19]. Results of computations with 
a large number of test problems have been reported in the 
cited papers and in each case these compact schemes are 
found to yield highly accurate numerical solutions. 

Moreover, the accuracy improves rapidly, consistent with the 
local truncation error, as the mesh is refined. 
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The Navier- Stokes equations representing the two- 
dimensional steady flow of an incompressible viscous fluid 
are given in streamf unct ion- vorticity form as follows: 


a j »///3x 2 + 

3 2 ^/3y 2 = -$ 

(5) 

9 2 $/3x 2 + 

3 2 S/9y 2 - Re ( u 3$/3x + v 3$/3y) = 0 

(6) 

u = 3i///3y, 

v = -3^/3x 

(7) 


Here 4> is the streamf unction, $ the vorticity; u,v are the 
velocities; Re is the nondimensional Reynolds number. 

The streamf unct ion equation (5) is a special case of 
eq.(l). The fourth order compact approximation for this 
equation may be obtained by putting u = \p , f = -$ and 
p(x,y)=0, q(x,y)=0 in eq.(3): 

4[i//,+ \p 2 + 4/,+ + ^ s + «//,+ \p t - 20^o 

= S» + $ 3 + $ « + 8$ 0 ]/2. (8) 

The vorticity equation (6) is also a special case of eq.(l) 
and the fourth order approximation in this case may be 
obtained by putting u=$, f=0 and p(x,y)=-Re u(x,y), 

q{x,y)=-Re v(x,y) in eq.(3),(4): 

8 

I c $ - c $ = 0, where (9) 

i=l i i 00 
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Cj = 4 - Reh/4 ( 4u„ +3u ! +u 2 -u 3 +u 4 ) 

+ (Reh) 2 /8 [4u 0 2 + u 0 (Uj— u 3 ) + v,(u,-u 4 )], 

c 2 = 4 - Reh/4(4v 0 +v 1 +3v 2 +v 3 -v4) 

+ (Reh) 2 /8[4v 0 2 + uJvj-Vj) + v 0 (v 2 -v 4 )], 

c 3 = 4 + Reh/4 ( 4u 0 -u i +u j +3u 3 +u 4 ) 

+(Reh) 2 /8[4u„ 2 - u 0 (u 1 -u 3 ) - v 0 (u 2 -u 4 ) ], 

c 4 = 4 + Reh/4 ( 4v 0 +v i -v 2 +v 3 +3v 4 ) (10) 

+ (Reh) 2 /8 [4v 0 2 - u 0 (v 1 -v 3 ) - v 0 (v 2 -v 4 )], 

c 5 = 1 - Reh/2(u 0 +v 0 ) - Reh/8 ( v 3 +u 2 -v 3 -u , ) + (Reh) 2 u 0 v 0 /4, 

c 6 = 1 + Reh/2 ( u 0 -v, ) + Reh/8 (v^Uj-Vj-u, ) - (Reh) 2 u 0 v 0 /4 , 

c 7 = 1 + Reh/2(u 0 +v 0 ) - Reh/8(v 1 +u 2 -v J -u 4 ) + (Reh ) 2 u 0 v 0 /4 , 

c, = 1 - Reh/2(u„-v 0 ) + Reh/8 (v 1 +u 2 -v 3 -u 4 ) - (Reh) 2 u 0 v„/4 , 

c„ = 20 + (Reh) 2 (u 0 2 + v 0 2 ) - Reh(u!-u 3 ) - Reh(v 2 -v 4 ). 

The velocities u,v at a grid point (x,y) are calculated 
from the discrete approximations of eq.(7). The typical 
second order central difference approximations for the 
velocities are: 

u„ = (^ 2 - ^«)/ 2h, v 0 = (ypi- ^i)/2h. (11) 

High order approximations for the velocities u, v can also 
be defined. We [10] earlier derived some high accuracy 
compact approximations for the gradients of the solution of 
Poisson equations. As the streamf unct ion equation (5) is a 
Poisson equation in \p , high accuracy approximations for the 
gradients 9\^/9x, 9i///9y can be obtained from [10] and 

corresponding approximations for the velocities are given 
below. (For details of the derivation, see the cited 
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reference.) The following approximations are compact and 
have a local truncation error of order h 4 . 

u 0 = (^i~^«)/3h + ( ^ 5 +v// 1 7 8 ) /12h + h($ 2 -$ 4 )/12 

( 12 ) 

v 0 = U,-iM/3h - (*,-*.-*,+*, )/12h + h(5,-5i)/12. 

MODEL PROBLEM 

As a model problem, we consider the steady flow of an 
incompressible viscous fluid in a square cavity (0<x<l, 
0<y<l). The flow is induced by the sliding motion of the top 
wall (y=l) from right to left, and is described by the 
Navier- Stokes equations (5)— (7) . The boundary conditions 
are those of no slip: on the stationary walls u=0 and v=0; 

on the sliding wall u=-l and v=0 (see Figure 2). 

A large number of investigators have used this model 
problem to test new schemes and solution methods, (see for 
example [4,5,6,7,8,9,11,17,18,20,22] and references given 
therein). Highly accurate benchmark solutions of this 
problem are available in the literature. In particular, Ghia 
et al. [7] obtained highly accurate solutions using 256x256 
grids for 100<Re<10000 . Schreiber and Keller [17] solved 
this problem using a continuation method on a sequence of 
grids including an 180x180 grid; Goodrich and Soh [8] used a 
streamf unct ion algorithm on a 65x65 grid. These solutions 
facilitate comparison and assessment of new solution 
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techniques. Experimental and numerical work on the three- 
dimensional cavity has been reported by Freitas et al.[4,5]. 
In order to solve the driven cavity problem, we replace 
the Navier- Stokes equations (5)-(6) by the finite 
difference approximations given in eq.(8)-(10). The 
velocities, defined in eq.(7), are calculated using either 
the second order approximations (11) or the fourth order 
compact approximations (12) in order to compare the 
effectiveness of these approximations. The unit square is 
covered by a grid of uniform mesh width h (h=l/N). The 
discrete approximations (8), (9) are written at each of the 
(N-l) 2 interior grid points. Zero values are prescribed for 
^ on the boundary; vorticity $ on the boundary is obtained 
using the Jensen formula [11,16]: on the stationary walls, 

we define $ 0 = (-8^ + ^ 2 )/2h 2 ; on the sliding wall y=l, we 
define $ 0 = (- 8^1 + \£ 2 )/2h 2 + 3/h. Here the subscript 'O' 

denotes a grid point on the boundary; the points '1', '2' 

lie inside the flow region such that the points 0,1,2 all 
lie on the straight line normal to the boundary (see Figure 
3). These approximations have local truncation error of 
second order. Higher order approximations could also be 
defined for obtaining boundary values of $; it is 

anticipated that the impact on the accuracy of the computed 
solutions would be marginal. These and other boundary 
approximations for vorticity were studied in detail by Gupta 
and Manohar [11] (see also Gajjar[6]). 
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An inner- outer iteration procedure is utilized to obtain 
the numerical solutions. At each outer iteration, the linear 
systems from the discrete streamf unct ion and vorticity 
equations are solved iteratively. We solved these linear 
systems using point- S.O.R. iteration with the relaxation 
parameters 1.7 and 1.2, respectively. The inner iterations 
are allowed no more than a preset number (usually 10) of 
iterations. As the fourth order approximations of the 
vorticity equation are stable, the inner iterations are 
convergent for all values of Re. 

We obtained numerical solutions of the driven cavity 
problem for Reynolds number ranging between 1 and 2000. The 
solutions were obtained on a 21x21 grid (h=l/20) and a 41x41 
grid (h=l/40). All iterations were started with zero initial 
data and were terminated when the maximum difference between 
successive approximations of \p , $ was smaller than 10' 4 . 
The computations were carried out on an IBM 4381 at The 
George Washington University and on a CRAY XMP24 at NASA 
Lewis Research Center. 

In Table 1, we present the representative parameters of 
the driven cavity problem for the 41x41 grid, obtained using 
the fourth order approximations (12) for the velocities. 
This table contains the values of <J/, $ at VC (VC= center of 
the primary vortex), the values of \f/ at the centers of the 
secondary vortices in the lower corners, and the value of $ 
at the mid-point (.5,1) of the moving wall. These 
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parameters are the major indicators of the accuracy of the 
computed solutions [11] and are quoted by most authors. 


Primary Vortex Secondary Vortices 


Re \p 5 



VC 

VC 

1 

.100027 

3.33906 

10 

.100029 

3.35029 

100 

.103463 

3.28572 

400 

.112814 

2.30247 

1000 

.111492 

2.02763 

2000 

.099586 

2.24579 

★ 

-.2091 (-5)= 

-.2091x10- 5 



* 

i — 1 

ID 

Xjn 

uvc 

DVC 


2091 ( -5 ) * 

- . 2100 ( -5 ) 

5.8637 

2212 ( -5 ) 

- . 2011 ( -5 ) 

5.8686 

1245 ( -4 ) 

- . 1747 ( -5) 

6.5505 

6512 ( -3 ) 

- . 1452 ( -4 ) 

10.0856 

1833 ( -2) 

- .1491 ( -3 ) 

15.9470 

2849 ( -4 ) 

- . 7597 ( -4 ) 

18.5790 


Representative Parameters of Driven Cavity 
(Fourth Order Velocities, 41x41 mesh) 

Table 1 


We also computed numerical solutions using the second 
order approximations (11) for the velocities. The solutions 
for Re=l,10 were found to be almost identical with those 
obtained using the fourth order approximation for the 
velocities. The representative parameters for second order 
velocities for Re>100 are given in Table 2. 



Pr imary 

Vortex 

Secondary Vortices 


Re 


$ 

* 

* 

$( .5,1) 


VC 

VC 

UVC 

DVC 


100 

.103263 

3.28369 

- . 1241 ( -4 ) 

- . 1742 ( -5 ) 

6.5641 

400 

.111151 

2.29561 

- . 7004 ( -3 ) 

- . 1367 ( -4 ) 

10.1538 

1000 

.107392 

2.01499 

- . 2108 ( -2 ) 

- .1384 ( -3 ) 

16.2462 

2000 

.088152 

2.37916 

- . 1368 ( -4 ) 

- . 7979 ( -4 ) 

20.6838 


Representative Parameters of Driven Cavity 
(Second Order Velocities, 41x41 mesh) 


Table 2 
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COMPARISON WITH EXISTING SOLUTIONS 


We now present a comparison of our solutions with the high 
accuracy solutions available in the literature (see, for 
example [7,8,17]). Qualitatively, our solutions exhibit the 
well known features of the driven cavity, including the main 
vortex in the central part of the cavity and secondary 
vortices in the lower corners. Figures 4-7 present the 
streamf unct ion and vorticity contours for Re=l and 100 using 
the 0(h 4 ) approximation for velocities. Figures 8-11 
contain the streamf unct ion and vorticity contours for Re=400 
using 0(h 2 ) and 0(h 4 ) approximations for velocities. 
Figures 12-15 contain similar contours for Re=1000. It is 
apparent that our streamf unct ion and vorticity contours are 
consistent with the published data [5,7,8,17,18]. 

Quantitatively, our solutions using the 41x41 mesh 
compare very well with the highly accurate benchmark 
solutions available in the literature. In Table 3, we 
compare the values of streamf unct ion ^ at VC and the 
location of VC with the results from [7,8,17] as applied to 
the driven cavity configuration shown in Figure 3. We note 
that the locations of the vortex center VC using the fourth 
order approximations for the velocities are within the 
cellwidth h (h=.025) of the reference data. 

Taking the results of Ghia et al.[7] as the benchmark 
solutions, we compute the relative errors of the values at 
VC for the solutions obtained by us, by Goodrich and Soh 
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Re 

Second Order 
Velocities 

Fourth Order 
Velocities 

Reference 

Data 


1 

.100027 
( .5, .775) 

.100027 
( .5, .775) 

.10006 [17] 

( .5, .7667) 

100 

.103263 

( . 375, .75) 

.103463 

( .375, .75) 

.103423 
( .3828, . 

[7] 

7344) 




.10330 

[17] 

400 

.111152 

( .45, .625) 

.112814 
( .45, .60) 

.113909 
( .4453, . 

[7] 

6055) 




.11198 

.11297 

[8] 

[17] 

1000 

.107392 

.111492 

.117929 

[7] 


(.475, .600) (.475, .575) (. 4687 ,. 5625 ) 

.11359 [8] 

.11603 [17] 


Value of \p at VC and location of VC 
(VC = Center of Primary Vortex) 

Table 3 

[8], and by Schreiber and Keller [17]. This data, given in 
Table 4, shows that our solutions obtained using the fourth 
order approximations for the velocities are either 
comparable in accuracy or are somewhat more accurate than 
the other results. 

The values of streamf unct ion </> at the center of the 
secondary vortex in the upstream corner UVC (in bottom left 
corner of the cavity, see Fig. 3) are given in Table 5. We 
also give the reference data from [7,8,17] for comparison. 
The values of this parameter obtained using the fourth order 
velocity approximations are much more accurate than those 
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Re 

Second Order 
Velocities 

Fourth Order 
Velocities 

Reference 

Data 

100 

.15% 

.04% 

.12% 

[17] 

400 

2.4% 

.96% 

.82% 

1.7% 

[17] 

[8] 

1000 

8.9% 

5.4% 

1.6% 

3.7% 

[17] 

[8] 


Relative Error in the Value of 4i at VC 


Table 4 


obtained using the second order velocity approximations; at 
large values of Re, the improvement becomes even more 
pronounced. 


Re 

Second Order 

Fourth Order 

Reference 


Velocities 

Velocities 

Data 


1 

- . 2091 ( -5 ) 

- . 2091 ( -5 ) 

-.247 (-5) 

[17] 

100 

- . 1241 ( -4 ) 

- . 1245 ( -4 ) 

- . 1254 ( -4 ) 

[7] 




- . 1320 ( -4 ) 

[17] 

400 

- . 7004 ( -3 ) 

- . 6512 ( -3 ) 

- . 6424 ( -3 ) 

[7] 




- . 5749 ( -3 ) 

[8] 




- . 6440 ( -3 ) 

[17] 

1000 

- . 2108 (-2 ) 

- . 1833 ( -2 ) 

- . 1751 ( -2 ) 

[7] 




- . 1892 ( -2 ) 

[8] 




- . 1700 ( -2 ) 

[17] 


Value of at UVC 

(UVC = Center of Upstream Corner Vortex) 
Table 5 


In Table 6, we give the values of $ at VC and compare 
with the available data. The agreement with the reference 
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data is quite good even though the location of VC 
substantially effects the values of this parameter; smaller 
values of h would locate VC more accurately and give even 
better agreement in the $ values at VC. In Table 7, we give 
the values of $ at the mid-point of the sliding wall y=l and 
the only comparison data available in [7], 


Re 

Second Order 
Velocities 

Fourth Order 
Velocities 

Reference 

Data 

1 

3.339 

3.339 

3.232 [17] 

100 

3.284 

3.286 

3.167 [7] 
3.182 [17] 

400 

2.296 

2.302 

2.295 [7] 
2.281 [17] 

1000 

2.015 

2.028 

2.050 [7] 
2.026 [17] 


Value of 

$ at VC 



Table 

6 


Re 

Second Order 
Velocities 

Fourth Order 
Velocities 

Reference 

Data 

1 

100 

400 

1000 

5.8637 

6.5641 

10.1538 

16.2462 

5.8637 

6.5505 

10.0856 

15.9470 

6.5745 [7] 
10.0545 [7] 
14.8901 [7] 


Value of $ at ( .5,1) 
Table 7 


In Table 8, we give the extreme values of the horizontal 
velocity u at the centreline x=0.5 near the bottom wall y=0. 
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Data from [7] and relative deviation of our results from 
this data is given for Re>100; also given is the relative 
error data from [8] for Re=400 and 1000. Similar data is 
given in Table 9 for the extreme values of the vertical 
velocity v at the centreline y=0.5 near the left wall x=0. 
We observe that our extreme u values (Table 8) are 
comparable to those of [8] although our extreme v values 
(Table 9) are somewhat more erroneous than those of [8]. 


Re 

Second Order 

Fourth Order 

Reference 


Veloc i t i es 

Velocities 

Data 


1 

.2065 

.2070 

- 


100 

. 2212 

.2223 

.2109 

[7] 


(4.9%) 

(5.4%) 



400 

. 3235 

.3288 

.3273 

[7] 


(1.2%) 

( .46%) 

(2.2%) 

[8] 

1000 

.3473 

. 3596 

.3829 

[7] 


(9.3%) 

(6.1%) 

(4.6%) 

[8] 


Extreme Value of Horizontal Velocity u 
at centerline x=0.5 near the bottom wall y=0 

Table 8 


In Table 10, we compare the values of and $ at VC 
obtained with the 21x21 and 41x41 grids and note that the 
values of these parameters rapidly approach towards the 
benchmark values when the grid is refined. The rate of 
convergence is somewhat slower when Re is large. It is 
expected that on further grid refinement, the solutions for 
large Reynolds numbers would also exhibit rapid convergence. 
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Re 

Second Order 

Fourth Order 

Reference 


Velocities 

Velocities 

Data 


1 

-.1678 

-.1685 

- 


100 

-.2263 

-.2289 

-.2453 

[7] 


(7.8%) 

(6.7%) 



400 

-.4022 

-.4203 

-.4499 

[7] 


(10.6%) 

(6.6%) 

(1.6%) 

[8] 

1000 

-.4188 

-.4522 

-.5155 

[7] 


(18.7%) 

(12.3%) 

(3.9%) 

[8] 


Extreme Value of Vertical Velocity v 
at centerline y=0.5 near the left wall x=0 

Table 9 



Re 

21x21 

solution 

41x41 
solut ion 

Reference 

Data 


1 

.099994 

.100027 

.10006 [17] 


100 

.103168 

.103463 

.103423 [7] 

VC 

400 

.101073 

.112814 

.113909 [7] 


1000 

.084237 

.111492 

.117929 [7] 



1 

3.0294 

3.3391 

3.232 

[17] 

$ 

100 

3.1112 

3.2857 

3.1665 

[7] 

VC 

400 

2.3784 

2.3025 

2.2947 

[7] 


1000 

3.0194 

2.0276 

2.0497 

[7] 


Comparison of 21x21 

and 41x41 

Solutions 



<//, $ at the Center of the Primary Vortex 
(Fourth Order Velocities) 

Table 10 
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DISCUSSION 


We note that for moderate values of the Reynolds number Re, 
the numerical solutions obtained using our high order 
compact difference approximations are highly accurate and 
compare well with the benchmark solutions available in the 
literature. This fact is remarkable also due to the fact 
that our solutions are obtained with a relatively coarse 
grid (h=l/40) whereas the benchmark solutions have been 
obtained with fine grids (with h as small as h=l/256). 

In Table 11, we give the number of iterations needed to 
converge to the required tolerance. As expected, for larger 
values of Re the convergence is slower; however, the 
convergence is faster with the fourth order velocity 
approximations than with the second order velocity 
approximations. We also give the CPU execution times for 
Re=1000 , 2000 on a CRAY XMP24 in this table. With Re=1000, 
the convergence with fourth order velocity approximations 
required almost 13% less CPU time than with the second order 
velocity approximations; with Re=2000 the discrepancy 
increased even further. Thus, the numerical computations 
using the second order velocity approximations have slower 
convergence and produce less accurate results than the 
computations using the fourth order velocity approximations. 

We conclude that the fourth order approximations for the 
Navier- Stokes equations do provide highly accurate 
solutions when coupled with appropriate high order 
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Re 


Second Order 
Velocities 


Fourth Order 
Velocities 


1 

157 


100 

353 


400 

516 


1000 

1248 

( 148 

2000 

>6200 

(>645 


157 

353 

509 

s) 1040 (127 s) 

s) 4266 (482 s) 


Number of Outer Iterations, 

CPU execution times 
to converge to 10" 4 
(CRAY XMP 24) 

Table 11 

approximations for the velocities. The rate of convergence 
of the outer iterations slows down considerably when 
computing for larger values of Re. This is an on-going 
difficulty with all solution methods. We are currently 
investigating alternative methods such as the multigrid and 
multilevel techniques [7,20,22] to obtain high accuracy 
driven cavity solutions for much larger values of Re, and to 
solve other problems of viscous fluid flow. The results of 
these investigations shall be reported in the future. 
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